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Abstract 

We study spectral and scattering properties of a spinless quantum 
particle confined to an infinite planar layer with hard walls containing 
a finite number of point perturbations. A solvable character of the 
model follows from the explicit form of the Hamiltonian resolvent ob- 
tained by means of Krein's formula. We prove the existence of bound 
states, demonstrate their properties, and find the on-shell scattering 
operator. Furthermore, we analyze the situation when the system is 
put into a homogeneous magnetic field perpendicular to the layer; in 
that case the point interactions generate eigenvalues of a finite multi- 
plicity in the gaps of the free Hamiltonian essential spectrum. 

1 Introduction 

The object of our interest in this paper is a spinless quantum particle living 
in a layer of a fixed width d with the Dirichiet boundary conditions and inter- 
acting with a finite number of point perturbations. An obvious motivation 
for this problem is to find a description for an electron in a semiconductor 
layer with impurities. However, such physical systems are in reality rather 
complicated objects which involve a crystal lattice with some alien atoms 
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and an electron gas, so one has to ask first whether such a simple model can 
reproduce the basic features known from experiments. 

It is well known that an electron in an ideally pure bulk semiconductor 
material can be modeled as a free particle with an effective mass m* which 
characterizes the relation between the energy and quasi- momentum at the 
Fermi level. Properties of the crystalline structure are thus expressed through 
a single material constant, which may be very different from the "bare" mass 
- recall that for GaAs we have m* = 0.067m e . 

There are two other assumptions in the "free" part of the model. The 
first is its one-particle character which neglects the interactions between the 
electrons. There are situations where the repulsion plays an important role, 
such as the Coulomb blockade in quantum wires. On the other hand, the one- 
electron model is known to work when the electron-gas density is sufficiently 
low. Another assumption is the neglection of spin which is also not entirely 
trivial; recall that spin-dependent effects in nanostructures have been studied 
recently - see [|l[] and references therein. In most situations, however, spinless 
electrons are a reasonable approximation. 

The next question concerns the way in which we model the impurities. 
Using again a certain idealization, we describe them by point interactions. 
This method proved rather useful in the last two decades and gave rise to nu- 
merous solvable models; our aim here is to add one more class to this family. 
Intuitively point interactions are understood as sharply localized potentials, 
but it is known that a sophisticated coupling constant renormalization is 
required to give this concept meaning in terms of a limit of scaled poten- 
tials H, Sees. LI, 1.5]. Mathematically such operators can be handled since 
they differ from the free Hamiltonian Hq just by a change of the boundary 
conditions at the interaction sites. However, the counterintuitive features 
of three-dimensional point interactions are reflected both in the slow way in 
which they found their place in the theory and in the fact that the parameters 
appearing in these conditions cannot be interpreted as 5 potential coupling 
constants but rather as the inverse scattering lengths corresponding to the 
point "obstacles" - cf. [Q, Chap. 1.3]. 

Scores of papers discussing point-interaction models in the Euclidean 
space, both for particles otherwise free or with a background regular po- 
tential, are summarized in the monograph ||. Only in the last decade the 
attention shifted to systems with point interactions restricted to a certain 
region of configuration space; the reason clearly was a wide collection of new 
physical phenomena observed in such spatially restricted systems, mostly 
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mesoscopic objects, but also electromagnetic waveguides, photonic crystals, 
etc. - see ||. Here, too, point-interaction Hamiltonians proved as a useful 
tool and yielded some unexpected results such as the existence of a chaotic 
behaviour in systems whose classical counterparts are integrable 

Today there are many papers treating point interaction in restricted areas; 
a bibliography is given in the introduction of |J. They typically put emphasis 
on the description of a specific model rather than a proper handling of the 
point interaction. Among few existing rigorous treatments of the problem it is 
the paper || which motivates the present study analyzing point interactions 
in an infinite planar strip with Dirichlet boundary conditions, together with 
similar systems. There are two ways in which the results can be generalized to 
dimension three. One is a straight Dirichlet tube in R 3 with a fixed compact 
cross section discussed in 0; it is a straightforward extension, apart of a 
different way of computing the regularized Green's function. 

In the present paper we are going to study a less trivial generalization, 
with point interactions situated in a flat layer with a Dirichlet boundary. The 
free system allows here again a separation of variables, so the free resolvent 
kernel and all the quantities derived from it such as eigenfunctions, etc., can 
be written by means of explicitly given series (in this sense models considered 
here are little "less solvable" than those in the full space when such quantities 
can be written in terms of elementary or special functions). 

Although the model description is simple, it covers many different situa- 
tions. For the sake of brevity we restrict ourselves in this paper to systems 
with a finite number of point perturbations in absence of a background po- 
tential leaving other cases to a sequel. We make an exception, however, by 
devoting a separate section to the case when the particle is under influence 
of a homogeneous magnetic field perpendicular to the layer. The spectrum 
of the unperturbed system is then changed completely consisting of infinitely 
degenerate eigenvalues which are sums of the Landau levels and the trans- 
verse eigenvalues; for "rational" combinations of parameters different Landau 
levels may lead to the same eigenvalue. A finite number of point perturba- 
tions then gives rise to a nontrivial discrete spectrum. 

Let us describe briefly the contents. The next section is devoted to the 
case of a single perturbation. We start from the definition of the point- 
interaction Hamiltonian by means of boundary conditions coupling gener- 
alized boundary values. After that we use Krein's formula to derive the 
explicit expression for the resolvent; it involves the regularized Green's func- 
tion which is given by a specific series as mentioned above - see ( 2.1 4j) . In 
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Section 2.4 we use this result to analyze spectral properties of such Hamil- 
tonians. The bound state energies are given by the implicit equation (|2.19|) , 
and it is just the limits of strong and weak coupling where we are able to 
write the explicit expressions for the leading term of the asymptotics. In 
both the extreme cases the eigenvalue behavior can be easily understood: in 
the strong-coupling situation it goes to — oo in the same way as if there were 
no boundaries because the corresponding eigenfunction is strongly localized, 
while in the weak-coupling case the eigenvalue approaches the threshold of 
the essential spectrum and the wavefunction is dominated transversally by 
the lowest mode. We also find that the eigenvalue decreases with the dis- 
tance from the layer boundary. In the last part of Section 2 we shall discuss 
the scattering in the presence of a perturbation. If there is a single point 
interaction, we can employ the symmetry of the problem with respect to ro- 
tation around the axis passing through the perturbation and perpendicular to 
the layer. The partial wave decomposition in the "longitudinal" coordinates 
shows that the only nontrivial contribution to the scattering comes from the 
s-wave, i.e. from states with the orbital momentum m — 0. Within this sub- 
space, the scattering problem is reduced to transitions between transverse 
modes; the final S-matrix then describes a coupling of the "open" channels, 
i.e. the transverse modes with the energies lower than that of the incoming 
spherical wave. We also derive the on-shell scattering operator which maps 
the incoming wave vector and transverse mode into the outgoing ones; the 
advantage of this approach is that it does employ the symmetry and allows 
for a generalization to the case with multiple perturbations. 

Section 3 extends the described analysis to any finite number iV of point 
perturbations. The technique remains the same, and since the difference 
between the two resolvents is of rank N, the essential and absolutely con- 
tinuous spectra are again preserved. On the other hand, the analysis of the 
discrete spectrum becomes more complicated. There are n eigenvalues, where 
1 < n < N, which are found by solving the implicit equation detA(2i) = 
with the N x N matrix A given by ( |3.8| ). The number n depends on the 
coupling strength. In the strong coupling limit there exist exactly A" eigen- 
values having the same asymptotic behavior as in the one-center case. On 
the other hand, for weak coupling we find only one eigenvalue approaching 
the threshold of the essential spectrum; in this sense our system exhibits a 
behavior typical for all weakly coupled Schrodinger operators. 

A new feature for systems with A^ > 2 is that they can posses eigenvalues 
embedded in the essential spectrum. This is possible, e.g., when the point 
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perturbations are placed symmetrically with respect to the layer axis and 
have the same (sufficiently strong) coupling: the corresponding eigenfunction 
cannot then contain contributions from transverse modes with the energy 
equal or smaller than this eigenvalue. We will show that this is true for 
embedded eigenvalues generally: their eigenvectors have to be orthogonal 
to the subspace spanned by the "lower" transverse modes. For N > 2 the 
system exhibits no longer a rotational symmetry, hence we cannot employ the 
partial-wave decomposition to describe the scattering. However, the second 
approach mentioned above is applicable here and we can derive again the on- 
shell scattering operator. It is similar to that of the one-centre case differing 
just by replacement of a single regularized Green's function by a sum of the 
elements of the matrix A - see ( |3.34| ). 

Section 4 deals with the situation when the layer is placed into a homo- 
geneous magnetic field perpendicular to its boundary. The Krein's formula 
is applicable but the free resolvent is substantially different from the non- 
magnetic case; this is reflected in the form of the essential spectrum which 
now consist the "sum" of the Landau levels and the transverse mode energies; 
it, of course, is preserved by a finite number of point perturbations. If there 
is a single perturbation we get exactly one eigenvalue in each spectral gap, 
i.e. between any two neighboring levels. In the strong and weak coupling 
limits it approaches the upper and lower endpoint of corresponding "free" 
gap, respectively. Only for the lowest gap we find a different behavior in the 
strong-coupling limit case; the eigenvalue goes to — oo with the same asymp- 
totics as in the non-magnetic case. Finally we present a generalization to the 
case of N point interactions analogous to the considerations of Section 3. 

2 A Single Perturbation 

2.1 The free system 

Consider an infinite layer E := M 2 x [0, d] with the coordinates denoted as 
x = (x,y) , where x = (xi,xz) G M 2 and y G [0,d]. We consider a single 
spinless nonrelativistic particle confined to E. Since the actual values of 
physical constants are not essential throughout the paper, we put h = 2m = 1 
and suppose that the free motion of the particle is governed by the Dirichlet 
Laplacian — A^. 

Recall that this operator can be defined for rather general domains in M. n 
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as the Friedrichs extension of an appropriate quadratic form |7], Sec. XIII. 15]. 
However, since the boundary of E is consists of two disjoint planes and has 
therefore the segment property, the operator acts simply as Hoip = —Atp 
on the domain of all ip of the local Sobolev space W 2,2 (T,) which satisfy the 
boundary conditions 

if)(x 1 0)=ip(x,d)=0 (2.1) 

for all xeM 2 - see again 0, Sec. XIII. 15]. 

We will make use of the fact that the "longitudinal" and "transverse" 
variables decouple in the free system. The state Hilbert space of our problem 
can be then decomposed into traverse modes, I/ 2 (E) = L 2 (R 2 ) ® 

because the functions Xn(y) '■— \J^ s ^ n (^') form an orthonormal basis in 
L 2 ([0,d]). The free Hamiltonian can be correspondingly written in the form 
of a direct sum, 

Jtv d 2 d 2 / TTTl \ ^ 

H„ = 0A.®I„, fc.^- M - s?+ ( T ). (2.2) 

Since the resolvent (H Q — z)' 1 of the two-dimensional Laplacian is known 
explicitly, the above decomposition yields in turn the free resolvent kernel 

oo 

G (x,y;x',y';z) = - ^ H { 1] (k n \x - x'\) Xn{y) Xniv'), (2.3) 



n=l 



where k n = k n (z) := y/z - (?f) 2 . 



2.2 Definition of a point interaction 

Our first goal is to construct a one-center perturbation supported by a point 
a := (a, b) with d£l 2 and b G (0, d). This can be done in a standard way 
Q. We restrict — A^, to functions which vanish in a neighborhood of a; the 
operator obtained in this way is symmetric but not self-adjoint and we look 
for the perturbed one among its self-adjoint extensions. Since the restriction 
has deficiency indices (1,1), the family of extensions is by the standard von 
Neumann theory |7|, Sec. X.l] characterized by a single parameter. 

What is equally important is that the perturbation is local, and therefore 
we can characterize the extensions by the usual boundary condition derived 
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in Chap. 1.1] for point interactions in R 3 . We introduce the generalized 
boundary values, 

L (ip,a) := lim \x — a\ip(x), Li(tp,a) := lim ip(x) , 

|i-3|-»0 |x-a|^oL \x — a\ J 

and require 

Lxty, a) -4iraL (ip,a) = 0. (2.4) 
For a fixed «6R this leads to the self-adjoint operator H(a, a) acting as 

(H(a, a)ip)(x) = -(Aip)(x) (2.5) 

for x a on the domain 

D(H(a,a)) := lip : — A-0 G L 2 (S \ {a}) and (2.1), (0) are satisfied} 

(2.6) 

in L 2 (S), where — A^ is, of course, understood in the sense of distributions. 
The family of self-adjoint extension includes also case which is formally given 
by a = oo, which means L (ip, a) = 0. It is easy to see that the corresponding 
H(a,a) is nothing else than the free Hamiltonian Hq. 



2.3 The resolvent 

As usual the spectral properties of the operator H(a,a) can be studied using 
its resolvent. Since H(a,a) and H have a common restriction with deficiency 
indices (1, 1), the kernel of the full resolvent can be obtained by means of the 
Krein's formula ||, Appendix A] 

,-rr, -s — * \ n -> \ i G (xi, a; z) G (a, x 2 ; z) 

(H(a,a)-z) {xi,x 2 ) = G {x 1 ,x 2 ; z) H — , (2.7) 

a — 4 (a; z) 

where 

n=l 1 1 

(2-8) 

is the regularized Green's function; we have employed here the fact that the 
resolvent singularity is the same as for the kernel of —A in L 2 (M 3 ) - see || 
Sec. 13.5]. The form of the denominator in expression (|2.7| ) follows from the 
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boundary condition (2^) applied to ip = (H( 



a, d) — z) 1 ip, where ip is an 
arbitrary vector from L 2 (£). However, the above definition does not give a 
practical way to compute £ (a; z) . To this end we use first K Q (z) = fH^\iz) 
and introduce K n := a/ (^f ) 2 ~ z = —ik n , then we have 



71=1 



(2.9) 



where g := \x — a\ and we have already put y = b. We use the asymptotic 
behavior /t„ « y as n ^ oo to write the identity 



d 



K {K n g) = K [n—) + K (K n g) - K [n— 



7TQ 



d 



(2.10) 



and to divide the function £ into two parts, £(a; z) = £i + £2; where 



- oo 

71=1 

1 oo 

6 == sL£h( 



7rn^ 



sin' 



7rn& 



V d 



'ring 
IT 



Ting 
~d~ 



cos 



2imb 



Aug 



we have used 2 sin 2 a = 1 — cos 2a. To deal with the first term we employ 
the asymptotic behavior of the Macdonald function ||, 9.6.13] which yields 



K (K n g) - K 



nng 
~d~ 



-In a/1 - z( 



d \2 



— (l + 0(^)) 
\7rn/ 



as g — > 0. It shows that the sum converges uniformly w.r.t. £> and the limit 
can be interchanged with the series giving thus 



6 



7TG? 



71=1 



irnb 



(2.11) 



The second part can be computed by means of the formula |10|, II, 5.9.1.4.], 

oo 



nx cos no 



n=l 



7T x - 



2V^ 2 + o 2 2 
1 



1 / , x 



^ L^/(2n7r - a) 2 + x 2 2mr 



2 ^V( 2wr + a) 2 + 



1 1 



x 2 2nir 



so introducing (3 :— K H '■= f and performing the limit /i — * , we get 
^ 2 = A^di~2p ~ ^ n(n 2 - (3 2 ) 1 ' 

n=l v ' 

The last sum equals using |10], I, 5.1.15.2.] expressed by means of the 
digamma function as ^(1) — |(^(l+/3) + ?/>(i— /?)), and since ^(1) = —7, where 
7 = 0.577... is the Euler number, and ip(l—/3) = ip (j3) + n cot(n (3), ip(l+/3) = 
ip{(3) + g, we arrive at 

Putting the results together, we get the sought formula 



«**) = -sS ta V 1 -K^)"* , ( ! r) + S5[' ,r+ *(3) + i 00t (T 

n=l 

(2-14) 

expressing the regularized Green's function in the form of a series. It is 
certainly more complicated than an expression of the corresponding quantity 
for the whole space in terms of elementary functions || Chap. 1.1], but it 
allows us to derive the needed properties of the function £ to compute the 
values of £(a; z) numerically. 

Remark 2.1 Notice the scaling behavior with respect to the change of the 
layer thickness, i.e. the formulae relating properties of the family S CT = 
M? x [0, da], o > 0. Here the dimension of the configuration space is decisive. 
While for two-dimensional system the scaling amounts to logarithmic shift 
of the function £ as shown in 0, in three dimension the transformation is 
multiplicative. We find easily that the situation is the same as for straight 
tubes in M 3 studied in [|], i.e. we have £(0^"; za~ 2 ) = a _1 £(a; z), where a? := 
(acr, ba) . This means, in particular, that the singularities of the resolvent 
which we will discuss below using equation (|2.19|) are related as follows: 

e a (a a ,a°) = a- 2 e(a,a) , a a := a~ l a . (2.15) 

Without loss of generality we put therefore d = tt in the rest of this and the 
next section. 
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2.4 Spectral properties 

The explicit form of the resolvent ( |2.7| ) allows us to derive information about 
the spectrum. Since its difference from (Hq — z)~ l is a rank one operator, 
the essential spectrum remains by Weyl's theorem J7|, Thm. XIII. 14] the 
same as for the free Hamiltonian H , i.e., we have a ess (H(a, a)) = [1, oo). 
At the same time, also the absolutely continuous spectrum is preserved, 
a ac (H(a, a)) = [1, oo), this time by Birman-Kuroda theorem |7|, Thm. XI. 9]. 

Next we would like to prove the absence of the singularly continuous 
spectrum. To this aim we have to check that the expression (if), (E z — E c )if)) 
for all if) G L 2 (E) is an absolutely continuous function for z from any interval 
(c, t) C (l,oo), in other words, that it can be written as an integral of a 
locally integrable function. Since a sc cannot be supported by discrete points, 
we may choose the interval (c, t) in such a way that it contains none of the 
points m 2 , m G N. For the spectral projection E t — E c to the interval (c, t) 
we employ the Stone's formula, 

1 /•* 

[if),{E t - E c )if)) = - Urn / (if),ImG(u + ie)if)) du; (2.16) 

we have used here the fact which we will establish a little later, namely that 
H(a,a) has no eigenvalues in (l,oo), and therefore the spectral projections 
to (c, t) and [c, t) are the same. The Green's function ( |2.7|) is analytic for 
z with Rez G (c, t) and ±lmz G (0,e), and furthermore, its limits when z 
approaches the real axis from above and from below exist and are continuous 
functions of z. Recall that by assumptions no thresholds are contained in 
(c, t) and £ has nonzero imaginary part for z ^ m 2 , 

Im£(a; z ) = — V s in 2 (n6) > , (2.17) 

n=l 



where we have used yfz — m 2 = i\/m 2 — z; hence the denominator of the sec- 
ond term in ( |2.7| ) cannot be singular. Consequently, the integrated function 
is bounded in (c, t) x [0,e) and by the dominated convergence theorem the 
limits can be interchanged with the integral giving 

i y* 

(if), (E t - E c )if)) = - / (if), Im G(u)if)) du. (2. 18) 

Jc 

The function under the integral is again continuous in the interval (c, t), hence 
it is integrable and the statement is proved. If t = to would be an isolated 
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eigenvalue embedded in the continuous spectrum and ip the corresponding 
eigenfunction, the above relation remains valid for t G (to—T],t ) U (t ,t Q +r]) 
with some rj > 0, while at the point to the l.h.s. should have a jump, which 
is clearly impossible due to the continuity of the integrated function. 

To determine the discrete spectrum, we have to find the poles of the 
resolvent. Recall that a perturbation which can be reduced to a self-adjoint 
extension of a common symmetric restriction with deficiency indices (1,1) 
can give rise to at most one simple eigenvalue in each gap of the spectrum 



TT] , Sec. 8.3, Cor. 1]. In our case it means one simple eigenvalue in the 
interval (— oo, 1). In view of the relation ( |2.7| ) one can find it by solving the 
implicit equation 

i{a-z) = a (2.19) 

for z G M. The series contained in the formula fl2.14|) converges for all 
z G K \ {n 2 : n G N}, because its terms decay like n~ 2 as n — > oo as we can 
see using the Taylor expansion of ln(l — £) to the first order. The remaining 
term £2 is independent of z and finite for any b G (0, tt). 

The value of £(a; z) is real for any z G (—00, 1). In particular, it is easy 
to compute 



£(a; z = Q) = 
Differentiating £(a, •) we get 



4vr 



7 + ip 



71 



+ ^cot(6) 



dz 2tt 



1 00 1 

— } — sin 2 (n&) > 



n=l 



(2.20) 



so the function is monotonously increasing for z < 1. Moreover, it diverges 
at both endpoints. At the continuum threshold we have 



f (a; z) 



71^ 



In y/l-z sin 2 (6) + 0(1) as 



(2.21) 



while on the opposite side we may employ a simple estimate 
1 



71- 



ln \/l — z sin 2 (6) + £ 2 —> —00 



as 



-00 . 



We will need a more precise asymptotics at large negative energies. Below 
we shall prove that 



£(a; z) 



4:71 



O e 



as 



z — > —00 



^2.22) 
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for any c < dist(a, bdfi) = ~ — |&— 1 1. In other words, the leading term 
corresponds to the analogous function for the point-interaction Laplacian in 
M 3 computed in 0, Chap. 1.1]; it corresponds to the heuristic concept that 
strongly bound states are well localized, and therefore not much influenced by 
the presence of the boundary. Another property of £(a, z) is its monotonicity 
across the halflayer, 

f(a; z) > z) if |6 - tt/2| < \b' - tt/2|, a = a'. (2.23) 

To prove it we employ the relation ( |2.8|) which yields 

oo 

t(a-z)-t{a!-z)= lim ^- V H^\k n {z)\x - a|)[sin 2 (n6) - sm 2 {nb')} • 

\x—a\— s-0 27T ^— ' 

Since sin 2 (n6) — sin 2 (ra6 / ) = sin(n(6 + 6')) sin(n(6 — b')) holds for < b' < b < 
7r/2 we arrive at 

f (a; 2) - ^(a; 2) = G (a, b + b';a,b- b'; z) ; (2.24) 

the monotonicity then follows from the positivity of free resolvent kernel - 
cf. 0, Appendix to Sec. XIII. 12]. This behavior is illustrated in Fig. |I[ 

This confirms the mentioned general conclusion: it follows from the stated 
properties of £(a; •) that the equation (|2.19|) has for any a G K a unique 
eigenvalue e(a,a) in (—00, 1] and that the function e(-,a) is monotonously 
increasing, 

e(a, a) > e(a', a) if a > a'. 

Furthermore, £(a, •) is a real-analytic function because for a fixed z\ < — 1 
it is expressed on a complex neighborhood of (— 00, zi) through a uniformly 
convergent series whose terms are analytic. It follows from the implicit- 
function theorem that e(-,a) is a C°° function - see [fl2| , Chap. XIV]. The 
function e(-, a) is also monotonous with 6 

e(-,a) < e(-,a!) if |6 - ir/2\ < \b' - ir/2\ . (2.25) 

The behavior of the eigenvalue is shown in Fig. 0. 

We are also interested in the asymptotic behavior of the eigenvalue in the 
limits of weak and strong coupling. In the former case we have 

( 2ir 2 a \ 

e(a, a) ~ 1 — exp j— as a — > 00 , (2.26) 

V sin {b)J 
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t) = 7t/6 

b> = Ji:/3 

b=7r:/2 
(z) 1/2 /(4tt:) 



-o.it 



d,:oc ^ 




0.20 



Figure 1: The function £(a; •) for three different positions of the point inter- 
action. The dash-dotted line is the leading term of the asymptotics (|2.22|) . 



where the symbol ~ means that for any e > and a large enough the 
eigenvalue can be squeezed between a pair of expressions from the r.h.s. in 
which 27r 2 is replaced by 2tt 2 ± e. On the other hand, the strong coupling 
asymptotics can be proved directly. By Dirichlet bracketing |7], Sec. XIII. 15] 
e(a,a) is for a < bounded from below by — (— Ana) 2 , and from above by 
the ground state of the Dirichlet Laplacian in a ball of radius c with the point 
interaction in the center. The latter is easily found: writing e(a, a) 
one has to solve the equation 

(-Ana) 2 = k 2 (1 + sinh™ 2 kc) . 

It yields the sought asymptotic behavior 

£ (a, a) = -16tt 2 « 2 + O (e ac ) as 



-K 



a 



-oo . 



(2.27) 



which justifies in view of ( |2.19|) a posteriori the relation ( |2.22|) . 

The formula ( |2.7| ) provides us with the (non-normalized) wavefunction of 
the bound state through the residue at the pole, ip(x; a, a) = G (x, a; e(a, a)), 
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b=7u/6 

b = 7T/3 
t>=7T/2 



Figure 2: The dependence of the eigenvalues e(a, a) on the parameter a for 
three positions of the point interaction. 

so we have 

oo 

ip(x; a, a) = — H^p yy/e{a, a) — n 2 \x — a\j sin(nfe) sin(ny). (2.28) 



n=l 



For a — > — oo we can write Hq(u) 
wavefunction is well localized, 



2_ e i(n-^/4) go we gee that 



ift(x; a, a) fa 



1 



3 — \/ 16n 2 a' 2 +n 2 \x—a\ 



n=l 



7r 3 \/l67r 2 a; 2 + n 2 \x — a\ 



sin(nfr) sin(ny) ; 

(2.29) 

this is illustrated in Fig. |3]. On the other hand, in the limit a —>■ oo the 
wavefunction for \x — a\ from a compact set behaves as 



, . _ . sin y 1 _ i . , . 

■?/>(£; a, a) « a — — m x — a smosmy 

sin o 7r 



(2.30) 



+ — if^ (Vl — ^ 2 1 a; — o| J sin(n6) sin(ny) . 



n=2 
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The wavefunction is dominated by the first transverse mode as Fig. |] shows. 

Let us finally return to embedded eigenvalues. We have excluded their 
existence away of the thresholds. If z = m 2 we use the equation Q2.7| ) again 
to check that the singularities at the r.h.s. cancel. In the vicinity of z = m 2 
the free resolvent kernel and the denominator a — £(a; z) behave as 

G (xi,x 2 ;z) = — ^\nVm 2 — z sin(rm/i) sin(rm/ 2 ) (l + O (V m 2 — z) V 
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a — £(a; z) = cH In \/ m 2 — z sin 2 (m&) 



7T 2 



where 



a := a — £ 2 H n > hi W 1 r sin 2 (n&) him sin 2 (m&) . 

7T 2 ^— ' V n 2 7T 2 

Then the full-resolvent kernel asymptotically behaves as 

( a -\ — -In \/m 2 — z sin 2 (mb) ] ^- In y/m 2 — z sin(myi) sin(m?/ 2 ) 

\ 7T 2 J TT 2 



X 



and cannot thus have a pole-type singularity at z — m 2 . 
Let us summarize the results obtained so far: 

Theorem 2.2 Let H(a, a) be defined by Qfi.dj) and ( \2.0j ) for d = tt; then 

(a) a ess (H(a, a)) = a ac (H(a, a)) = [1, oo) and a sc (H(a, a)) = 0. 

(b) For any a G R there is a single eigenvalue e(az,a) in (— oo, 1) which is 
increasing and infinitely differentiable w.r.t. a. The corresponding eigen- 
function is given by (\2.28j ). 

(c) The eigenvalue is by (\2.25 ) strictly monotonous across the halflayer, 



s(a, a) < e(a, a') if \b - n/2\ < \b' - ir/2\. 

(d) In the limit a — > oo the bound state behaves according to ( \2.26j) and 
( \2.3Q ). In the strong coupling case the eigenvalue asymptotics is given by 



( \2.2Tj ) and the eigenfunction is described by (\2.2£\) . 



(e) There are no eigenvalues in [l,oo). 



2.5 Scattering 

Since the Hamiltonian H(a, a) is invariant under rotations around the axis 
passing through the point a and perpendicular to E, we may simplify the 
treatment of stationary scattering using a partial-wave decomposition. We 
use the tensor-product representation 

L 2 (S) = L 2 ((0,oo) x [0,d];rdrdy) ® L^S 1 ) , (2.31) 
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where S 1 is the unit circle in IR 2 and r := \x\. This can be written as 



L 2 (S)= C/- 1 L 2 ((0 > oo)x[0,d|)®{y m }, 



^2.32) 



where U : L 2 ((0, oo)^x [0, d];rdrdy) — >■ L 2 ((0, oo) x [0, d]) is the unitary 
operator defined by (Uip)(r) := r 1 / 2 -?/;(r), the "spherical" functions l^(o;) : = 
(27r) _1 ' 2 e ime with = (cos 9, sin 0) form a basis in L 2 (S' 1 ) and the symbol 
{•} means as above the linear envelope. Shifting the point a to the origin of 
the polar coordinates by T a : (T a (f))(x,y) = (f>(x + a,y), we can decompose 
the corresponding free Hamiltonian H = — as 



[2.33) 



with the partial-wave operators 



2 



dr 2 



d 2 Am 2 
+ 



dy 



4t 2 



m G Z . 



(2.34) 



Their domains are given in a standard way J13], Sec. 5.7]; the only non-trivial 
part is the radial boundary condition at the origin. As in |2|, Chap. 1.5] none 
is imposed in "higher" partial waves, m^O, because the radial part of ( |2.34|) 
is then limit-point at zero |7], Appendix to Sec. X.l]. Consequently, this part 
of the Hamiltonian is trivial from the point of view of the point interaction. 
On the other hand, for m = we introduce the generalized boundary values, 

W){y) : = lim <j>(r,y)y/¥ , ^(y) := lim r^^r, b) - £ (cf>) {y)^ 1 ' 2 ] . 

r— »0 r— >0 

(2.35) 

The s-wave component of the free Hamiltonian is specified by the condition 
£o(<f))(y) = (for any y G (0, d)), while the s-wave component h$ of H(a, a) is 
given by the same differential expression ( |2.34| ) with the boundary condition 
at y — b changed to 



i 1 (4 > )(b)-4naW)(b)=0 



[2.36) 



It is also clear from this discussion that the eigenfunction of H(a, a) analyzed 
in the previous section exhibits a symmetry, ip(R((p)x; a, a) = ?/>(:r; a, a), 
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where R(ip) is the rotation of £ on an angle ip around the axis passing 
through the point a and perpendicular to E. 

Let us return to the scattering problem. As in 0, Chap. 1.5] the rotational 
symmetry means that the S-matrix part corresponding to partial waves with 
m G Z \ {0} is trivial, i.e. the unit operator. In distinction to || Chap. 1.5], 
however, the s-wave part (we denote it as S again to keep the notation simple) 
is still in general a complicated operator because the point interaction can 
couple different transverse modes. Its dimension depends on the number of 
the "open channels", i.e. of the transverse modes in which the particle of 
energy z can propagate; for d = ix the latter is [y/z\. Using the partial- wave 
operator ( |2.34| ) with m = it is easy to see that the function 

i/>W(£-k) = J {k n {z)r)xn(y) (2.37) 
1 v 

+ a _f/g. z \ ^ Yl H o\ k j( z ) r ) sin 0' 6 ) sin(n6)xi(y) 

with r = \x — a\ satisfies the boundary condition (2.4) and r l / 2 ipa)n{x] k) is 
a generalized eigenfunction of <g> ^— J^J with the eigenvalue z = k 2 . In 
the limit r — * oo it behaves as 



1>£l&k) « J—^—cos(k n (z)r-7T/A) X n(y) (2.38) 
y T^K n \z)r 

To find the matrix elements S n j one has to compare this with the asymptotics 
of the outgoing wave in the jth transverse mode expressed by means of the 
scattering phase shift. For the wave scattered back into the incident nth 
mode we have 



irk n (z)r 



COs{k n { Z )r - 7T/4) + — e i(km(z)r-v/4) sin 2( n fy 

a — £(a; z) 2n 



id 



n n 



cos(k n (z)r - 7r/4 + S nn (k)) 



irk n (z)r 
which yields 

7T a— £(a; z) 
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In a similar way scattering to the jth mode, n ^ j, requires the identification 

1_ e i(k 3 {z)r-*/i) sin / 6 x sin / - & x = e it nj (k) cos (k n ( z )r-n/4 + 5 nj (k)) ; 

a— 4(a; z) 27r 

together we get 

r. (k) r * sin(nfr) sin(j&) 

-S'nj(fc) = e 2 ^ (fc) = <y nj - + ; \ . 2.39 

We have to check that the obtained S-matrix is unitary, i.e. 

Y,S nj S sj = 5 ns . (2.40) 

3=1 

Using (|2.39| ) we can write 

^ S n jS s j = ^ .{^nj ~ 3nj)(S s j — 5j s ) + 2Re (S ns — 5 ns ) + 5 ns 
3=1 3=1 

7T |a-^(a;^)| 2 I 2vr 

so the desired property follows from ( |2.17 ). 

The scattering problem can be also described in another way - by means 
of a scattering operator in L 2 (S' 1 ) <8>L 2 ([0, d}). Applying (2771) to an arbitrary 
G £ 2 (£) we see that to any if) G D(H(a,a)) and a nonreal 2 there is 
if) z G D(Hq) such that 

^(5) = Vz(5) + G {x, a; ^)^(a). (2.41) 

If we choose, in particular, ip £ z (x) = e tk "( z ) UJX ~ e \ x \ 2 x n (y) with a; a unit vector 
in M 2 , then the corresponding ip e (x) G D(H(a,a)) satisfies the equation 

{{H(a,a)-z)if) £ )(x) = Ae[l-e\x\ 2 + ik n ujx] if) e z {x) . (2.42) 

The r.h.s. converges in the L 2 sense as z approaches the real line and the 
resulting if) 6 belongs to D(H(a,a)) for a fixed z G [l,oo). Furthermore, the 
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pointwise limit exists as e — » 0+ and equals 

i^aJx; k n (z)u) = e ik ^ x xn(y) (2.43) 



e ik n (z)uja ^ x 



+ — 7T J2 H o\ k j( z )\ x - a l) sin (i & ) s Hnb)xAy) ■ 

a - f (a; 2) 27r ^ 

The function defined by the r.h.s. is locally square integrable, solves the 
equation (H(a, a)—z)ip = and satisfies the appropriate boundary condition, 
i.e. it is a generalized eigenfunction of H(a, a). Let us expand it into partial 
waves. We know already the s-wave eigenfunction fl2.37|) , the remaining 
ones are trivial, ^^(x; k n (z)) = Jm(kr)x n (y) if m 7^ 0. Here again, this 
expression describes an eigenfunction in the Hilbert space fl2.31|) , 

multiplying it with r 1//2 one obtains an eigenfunction in the Hilbert space 
(|2.32[) . Using the known identity 

00 

e <M*-a) = 27r i m J m {k\x-a\)Yju)Y m {uJ x _ a ) (2.44) 

m=— 00 

where Y m are the functions introduced above and 0J x -a '■= ifEfr we get 
e-* 1 ^^^; K(z)uj) = ^S(x; k n (z)) 

Components of the on-shell scattering amplitude (f a (k n (z), ou x , oj))- are then 

given by Km.\ x \-+oo,x/\x\=u>* \x\^e~ tk ^ [ijj a!n (x; k n (z)u) - e tknLUX Xn(y)], specifi- 
cally its part corresponding to the j'th transverse mode is 

(f a (k n (z)^)) jn = sinQ-ft) sin(nft) 

and the on-shell scattering operator S a (z) on L 2 (S' 1 ) ®L 2 ([0, d\) has the form 

S a {z) = I+ l - V " D(n y" 5) ( e -*»W<>y X», •) e^ W( *%JG • (2-46) 

It follows from Q2.17| ) that the denominator is nonzero in [1, 00). However, we 
have argued above that £(a; •) can be continued analytically to the complex 
plane where zeros exist in general. In the weak coupling case there is one 
resonance pole of S a (z) close to the threshold of each higher transverse mode 
similarly as in the two-dimensional case |J. 
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3 A Finite Number of Point Interactions 



3.1 Boundary conditions 

Consider now a finite number N of point interactions and suppose that their 
positions are ctj = (aj,bj), where cij 6 M 2 and bj G (0,7r). For the sake 
of brevity we denote a = (ai, . . . , a at) and a = (cti, . . . , a^). The way to 
define a point interaction is the same as above; now we have N independent 
boundary conditions 

Li(-0,%) - 4irajL (ip,aj) = 0, j = 1, . . . , iV. (3.1) 



The Hamiltonian H(a,a) is given again by the formulae ( |2.5|) and Q2.6|) with 



the boundary condition ( p.4|) replaced by (|3.1|) and a understood in the sense 
mentioned above. Any of the point interactions can be switched off when 
corresponding coupling constant aj is formally put equal to infinity. 

3.2 The resolvent 

We again use the Krein formula to find the resolvent kernel. Since the defi- 
ciency indices of the operator obtained by restriction of the free Hamiltonian 
to the set of functions which vanish at the vicinity of the points a are equal 
to (N, N), the r.h.s. is a rank iV operator expressed in terms of vectors from 
the corresponding deficiency subspaces, 

(H(a,a) - zy 1 (xi,x 2 ) 

N 

= G (xi, x 2 ; z) + ^ ^jk{®, a; z) G {x x , af, z) G (a k , x 2 ] z) . (3.2) 
j,k=i 

Applying this to an arbitrary vector of L 2 (E) we get 

N 

if)(x) = t/) Q (x) + 2J A ifc (a, a; z) G (x, z) ^ (ofc) (3.3) 
j,k=i 

with ipQ G D(H ). The generalized boundary values are 
N X-k 

L (ip,a m ) = ^ ^-5 jm ip (a k ) , (3.4) 
j,k=i 71 
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N , 

Li(^,a m ) = Vo(a m ) + J XjkSjm lim ( G (x, a^z) 



47r I x — a, r 



iV 

+ ~ $jm)G (a m , a,-; z)^{a k ) • (3.5) 

j,k=i 

The limit contained in the expression of Za(/0, a m ) equals to £(a m ; 2). After 
substituting these boundary values into (|3.1|) we arrive at the conditions 

N 

ipo{a m ) + ^ ^ k l S jrn(£(a m ; z) - a m ) + (1 - 5 jm )G (a m , %; z)] ^o(afc) = , 
j,fc=i 

(3-6) 

which should be satisfied for an arbitrary vector i/jq belonging to D(Hq), i.e. 
for any A-tuple (^o(oi)) ■ • • > V'o(ajv))- This is possible only if the expressions 
in the square brackets are up to the sign elements of the matrix inverse to 
A(a, a; z), in other words, if the coefficients are 

A(a, a; z) = A (a, a; z)^ 1 , (3.7) 

where 

Aj-j := aij — £(aj',z) (3.8) 
= ^ + ^ f> \Z 1- ^ sin 2 (r%) _ -L T + ^ ^ + I cot(6 i 



Aj m :— —Go(aj,a m ;z) 



~ ^ #0^ (v^ - ™ 2 l a j - a ™l) sin(nfej-) sin(n6 m ) , j ^ m. 



n=l 



The Green's function value Go(aj,a m ;z) is finite for any pair of mutually 
different vectors Sj and a m . When the point interactions are arranged verti- 
cally, i.e. a,j = a m , the expression through Hankel's functions is useless and 
the corresponding nondiagonal element A jm can be alternatively written as 



A jm = — In Wl - 4r sin(nbj) sm(nb m ) (3.9) 

n=l 



To derive this expression we employ the argument analogous to that leading 



to the value of the function £(a; z) in Section 
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3.3 The discrete spectrum 

Since a finite rank operator is both compact and trace class, the argument 
presented at the opening of Section 2A remains valid. In other words, a finite 
number of point interaction changes neither the essential nor the absolutely 
continuous spectrum, a ess (H(a, a)) = a ac (H(a,a)) = cr ess (H ) = [1, oo). 
The singularly continuous spectrum is empty because the proof given in 
Sec. ^]4| can be used here again. 

The discrete spectrum is determined by poles of the resolvent, which 
occur if the coefficient matrix (Xjk)~ becomes singular. This leads to the 
condition 

det A(a,a;z) = 0. (3.10) 

To find the eigenfunctions we use the procedure from 0, Sec. II. 1]. Suppose 
that H := H(a, a) satisfies the equation Hp = zip for some z G K and pick 
an arbitrary z' G g{H) . Then in accordance with ( |3.3| ) there is a function 
■00 £ D(Hq) which makes it possible to write 

N 

p(x) = ipo(x) + ^2d j G (x,a j ; z') , (3.11) 

3=1 

with the coefficients dj := Yl!k=ii^ z ') 1 ) jk^^k) ■ We also see that 

(H - z')^ = (H- z')p =(z- z')ip. (3.12) 

Applying the resolvent (H — z')^ 1 to the last identity we arrive at the ex- 
pression 

N 

(h - z'y 1 ^ + ^2 d j( H o - z 'y lG o(-, z) 

3=1 

which allows us to find the action of (H — z) at the vector 

N 

{H - 2)Vo = (H - z')% -{z- z')i> = (z- z') d jG (; Sf, z') . (3.14) 

3=1 

If z < 1 , the resolvent (Hq — z)~ l exists and may be applied to the last 
relation, giving by means of the first resolvent identity 

N 

tp = ^d i [G (-,a J ; z) - G (-,a i ; z')\ . (3.15) 

3=1 



ip = (z- Z 



(3.13) 
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Substituting this into ( |3.11| ) we get an expression for the eigenf unction, 



N 

(p(x) = djGo(x, dj] z) , (3.16) 
3=1 

where it remains to determine the coefficients. The equation ( |3.15| ) taken at 
the point x = 3j can be rewritten using the components of matrix A as 

N 

ipoia-j) = ^d m (A(a, a; z') jm - A(a, a; z) jm ) . (3.17) 

m=l 

We already know that d m = '^2^ =1 (A(z l )~ 1 ) m k4'o(d k ); taking the inverse we 
get ipo{dj) = Y2 m =i ^{ z ')jmdm- In combination with (|3.17|) this yields 

N 

H^jmdm = 0, (3.18) 

m=l 

i.e., d := (di, . . . , d^) has to be an eigenvector of the matrix A(a, a; #) corre- 
sponding to zero eigenvalue. The corresponding system of linear equations is 
solvable under the condition ( |3.10| ), and the sought eigenf unctions of H(a, a) 
are given by the formula ( |3.11|) . 

It remains to check that the equation ( |3.10| ) can have a solution in 
(— oo, 1). Let us start with the limiting situations of strong and weak cou- 
pling. We know from (|2.22| ) how £(0,-; z) behaves as — -^p as z — > —00. At 



the same time, the nondiagonal part of matrix A vanishes in the limit in view 
of the asymptotics 

Hq (V z — n 2 \dj — a m \) = H^\iyjn 2 — z\aj — a m \) 



VK^ri 2 — z \ dj — a r 



e 



— sjri 2 — z\ dj— a m \—iir/4 



for z — > —00. This argument is not applicable if aj = a m for j 7^ m. 
Nevertheless, the nondiagonal matrix elements are up to the sign equal to 
Green's function values, and thus they are bounded. In this way we get 



A(a, a; z) = 1 + 0(1) as z -> -00 (3.19) 

4:71 
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with the coefficient azj included into the error term. On the other hand, 



A(a, a;z) = — \nVl - z M + 0(1) as z -> 1- , (3.20) 

TV 2 

where M := (sinfejSin& m ) . m=1 . This matrix has zero eigenvalue of mul- 
tiplicity N — 1 and the positive eigenvalue X)*Li s i n2 corresponding to 
the eigenvector (sin&i, . . . , sin&jv)- The latter is more important, because it 
means that one eigenvalue of A (a, a; z) tends to — oo as z — > 1— . Further- 
more, the eigenvalues of A (a, a; z) are continuous functions of z, so comparing 
the last claim with ( |3.19| ) we find that at least one of the eigenvalues crosses 
zero for some z, i.e. that H(a, a) has at least one isolated eigenvalue. 

We may ask whether there some of the eigenvalues may be degenerate. 
For the sake of brevity we rewrite the matrix A (a, a; z) as 

A(a, a; z) = 5 jm (aj - ^(z)) + (1 - 5 jm )g jm (z), (3.21) 

where £j(z) := £(aj]z) and gj m (z) := —G (aj,a m ;z). Since all gj m are 
negative, the maximum possible degeneracy is AT — 1 which means that for 
a pair of point interactions the discrete spectrum is always simple. Let us 
consider the case N = 3. Put z < 1, 01,3 = (±a, b\) and 0,2 = (0,62)- If 
b 2 — > then the value gi2 = (723 approaches zero being thus smaller than 
#13 for 6 2 small enough. On the contrary, if 6 2 = b% we obtain the opposite 
inequality; this follows from the expression (|2.3j ) and from the monotonicity 
of the Macdonald function K (u) with a positive argument - see || 9.6.24]. 
Hence there exists a 6 2 e (0, 61) such that all the non-diagonal elements of the 
matrix are the same, (713 = = #23- Choosing ctj which satisfy cxj — ^j = g\2 
we obtain a matrix of rank one, i.e. z is an eigenvalue of multiplicity two. 

3.4 Embedded eigenvalues 

We have shown in the previous chapter that a single point interaction cannot 
produce eigenvalues embedded in the continuous spectrum. This is no longer 
true if N > 2 as the following example shows. 

Example 3.1 Consider a pair of perturbations with the same a placed at 
d\ = (0,0,6) and a 2 = (0, 0, n — b). We can divide the eigenvalue problem 
into symmetric and antisymmetric parts with respect to the plane {(x, 7r/2) : 
x G M 2 }. We obtain properties of the antisymmetric part by scaling the one- 
center problem: substituting a — \ into the relation ( |2.15| ) we see that the 
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antisymmetric part has a single eigenvalue which tends to 4 as a 
it is embedded in the continuous spectrum for a large enough. 



oo, hence 



Thus we cannot exclude existence of embedded eigenvalues in general. 
We can, however, prove a weaker result. In the example the symmetry was 
essential, which means in particular that the eigenfunction is dominated by 
the second transverse mode. We will show that in general any eigenvalue 
z > 1 cannot contain contributions from transverse modes with n < [y/z\. 
Suppose that H(a,a)ip = zip for z > 1. We employ the relation fl3.11|) and 
take tpQ in the form i/)q(x) = Y^m=i 9n{x)Xn{y), where g n G L 2 (R 2 ). Substi- 
tuting this ipo into ( |3.11| ) and using the fact that {Xn} is an orthonormal 
basis in £ 2 ((0, 7r)) we get 



d 2 



d 2 



z + n 



9n(x) 



N 



3=1 



3 1 



for n — 1, 2, 



The Fourier-Plancherel operator transforms this into 



[|p| 2 - z + n 2 ] g n (p) 



\z - z' 
2tt 



N 



-ipaj 



3=1 



\p\ 2 — z' + n 2 



(3.22) 



If g n e L 2 



then g n should also belong to L 2 



It is not possible if 



z > n 2 and the r.h.s. of ( |3.22| ) is nonzero at p n co, where p n = V z — n 2 and 

1 has no 



z' + n 2 ) 



w is a unit vector in R. Recall that the factor 
singularity, because z' e p(H) by assumption. 

To avoid the singularity of g n at \p\ = p n we have to require 



N 



^djXnib^e 



ip n UJCLj 



(3.23) 



3=1 



for an arbitrary unit vector u) from M 2 . If all the a/s are different it follows 
that djXn{bj) = for each j. If some of them are the same the condition 
changes to J2j d jXn{bj) = where j runs through the a/s which coincide. In 
both cases g n is identically zero for n 2 < z. 
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Consider now an arbitrary g G L 2 (IR 2 ) and n 2 < z. Using and ( |2.3| ) 

we arrive at 



N 

(9Xn,<p) = (9, 9n) + ^dj Xn(bj)(g, -H^ 1] (k n (z')\x - aj\)) . (3.24) 

3=1 

The first scalar product is equal to zero because g n is zero. If all the a/s 
are different then djXn{bj) = for all j and the whole r.h.s. is equal to 
zero. If some a/s are the same, the condition ^2f =1 djXn(bj) = leads to 
(gXn, V 9 ) = again. The conclusion holds true for all n = 1, . . . , [-y/i] what 
we have set out to prove. 



3.5 The limits of strong and weak coupling 



Of the two extreme situations, consider first the strong coupling. One can 
write the matrix A(a, a; z) in the form 



A(z) 



An 



6 



/ + 



a,- + 



An 



-i 



S jk A(z) 



(3.25) 



where A(z) is the remainder matrix, which is independent of atj and has a 
bounded norm as z — > — oo. For an arbitrary finite interval J C (— oo, 1) 
we can always choose a/s large enough negative that no eigenvalues are 
contained in J and the two matrices in the above product are regular in 
J. It means that the roots of equation fl2.19 ) come from the region where 
A(z) is dominated by its diagonal part. Then there are exactly N eigenvalues 
(including a possible degeneracy), which behave asymptotically as 



€j(a, a) ~ — 16n 2 a 



as 



max a,- 

1<3<N 



— oo. 



(3.26) 



As we expect the eigenfunctions in the strong limit are strongly localized 
and only slightly influenced by the other perturbations. The eigenfunction 
localized at Oj has the same form as in the Section |2.4| where we put a — 3j. 

Consider now on the contrary that all the point interactions are weak, 
i.e. that all the <x,'s are large positive. Denoting A := diag(«i, . . . , a^) we 
can write 



A(z) = A + — In y/l-z M + A(z) 



7T Z 



(3.27) 
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where A(z) is a remainder matrix, which is again independent of a with its 
norm bounded, ||A(z)|| < C$ for z G (z ,l). In Section [T3] we found that 
A (z) was asymptotically a rank-one operator on C N . Hence only one isolated 
eigenvalue of H(a, a) exists in this asymptotic situation. 

To find the leading term of the asymptotic expansion, we have to solve the 
spectral problem for matrix M — rjA — r]A(z), where r] := — 27r 2 (ln(l — z))~ x . 
The largest eigenvalue of this matrix satisfies Hn(v) — A*jv(0) — (Ca+a + )rj, 
where a + := maxaj, while for j = 1, . . . , iV — 1 we have Hjijf) < (Cg— Q-)rj, 
where a_ := min<x,-. Since jUjv(O) > we see that for a_ large enough the 
condition has just one solution for rj > 0. One can check directly that without 
A (z) the condition is satisfied for 77 = £\ aj 1 sin 2 bj. Thus 77 = Oi^aZ 1 ) and 
the eigenvalue expansion follows. Hence the bound-state energy in weak- 
coupling case behaves as 

e(a,a) « 1 -exp |-2tt 2 J (3.28) 

as minjxjXTvOj — > 00. Since the eigenvector of matrix M corresponding to 
the nonzero eigenvalue is (sin & x , . . . , sin the asymptotic expression of the 
eigenfunction for x from a restricted part of S is 



~ sinw — — , ,, — — - > sm bj In x - 

^ v ' a I v-\iV sin 2 bj yj-2 1 

\ 2^j=i j=l 
H — - ^^sin(ra/) ^^sin&j sm(nbj) K {^Jn 2 — l\x — dj\j . (3.29) 

71 n=2 j=l 

Let us summarize the spectral properties of our Hamiltonian in the N- 
center case derived in the previous three paragraphs. 

Theorem 3.2 Let H(a, a) be defined by ( \2.3j ) and t \2.dj) for d = n, where 



a — (Si, ... , &n) and a = (a±, . . . , a at) and the boundary condition (\2.4) is 
replaced by ( \3. 1\) ; then 

(a) a ess (H(a,a)) = a ac (H(a,a)) = [1, 00) and cr sc (H(a,a)) = 0. 

(b) For any a G M. N there is at least one eigenvalue e(a, a) in (—00, 1). The 
maximum number of eigenvalues is N with the multiplicity taken into ac- 
count; the maximum multiplicity is N — 1 . The corresponding eigenfunction 
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is given by ( 3.16Q , where the coefficients d m , m = 1, . . . , N , are components 
of a vector solving the equation Ad = with the matrix A given by ( \3.8{ ). 

(c) In the limit mini<j<N ctj — > oo the bound state wave function behaves ac- 
cording to l \3. 2$ ) and ( \3. 2$ ). On the other hand, in the strong coupling case 
there are exactly N eigenvalues whose asymptotics is given by ^3.2b\ ); the 
corresponding eigenfunctions are strongly localized around the points x = dj 
and given by l \2. 29[ j with a = dj and a — atj . 

(d) If an eigenvalue z G [l,oo) exists, the corresponding eigenvector is or- 
thogonal to the subspace 0^ L 2 (R 2 ) <g> {Xn}- 



3.6 Scattering 

Comparing to the one-center case, the Hamiltonian H(a, a) with a finite 
number of perturbations loses in general the invariance with respect to rota- 
tions around an axis perpendicular to E. Hence we cannot employ here the 
partial wave decomposition and we turn directly to the "closed form" of the 
on-shell scattering amplitude and on-shell scattering operator. By (|3.3|) , to 
any ip G D(H(a, a)) and a nonreal z there exists ip z G D(H ) such that 

N 

ip{x) = ip z (x) + ^ ^jk{u, a; z) G (x, ay, z) ip z (a k ). (3.30) 
j,k=i 

We take again i/Jq(x) = ^M 2 )^-^! 2 -^^) f or ^ where u is a unit vector in 
IR 2 . Denoting the l.h.s. of ( |3.3U| ) as ip £ we have if) 6 G D(H(a,a)) and 

((H(a,d)-z)ij £ ) (x) = 4e[1-e\x\ 2 + ik n ux) ip E z (x) . (3.31) 

The r.h.s. converges in L 2 (M. 2 ) as z approaches the real line and the tf) £ 
belongs to D(H(a, a)). The pointwise limit e — > 0+ exists and equals 

N 
j,k=l 

(3.32) 

The limiting function is locally square integrable and it thus is a generalized 
eigenfunction of H(a, a). 

The components (f a {k n {z),u) x , u))) mn of the on-shell scattering amplitude 
are then given by the part of the following expression corresponding to the 
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outgoing mth transverse mode, 



|x|— >oo,x/\x\=u) x 



It yields 



Mr/4 • 



(f a (K(z),u; x ,u;)) mn = -^^e^A^oiz)^ 

x sin(mfej) sin(nfejfc) , (3.33) 
and the on-shell scattering operator S^z) on L 2 (S' 1 ) <8> L 2 ([0,rf]) is 

2- 2 



fc,jf=l m,n=l 

x e -ifc m ( 2 )(.)a, Xm _ (334) 

As in the one-center case, resonances are determined by the poles in the 
meromorphic continuation of the matrix-valued function (Xj k (a, a; •)). 



4 A Layer in Magnetic Field 
4.1 The free Hamiltonian 

In this section, the layer E = R 2 x (0, d) is placed into a homogeneous 
magnetic field B = (0,0, B). As usual the vector potential generating this 
field can be chosen in different ways, e.g. we can employ the symmetric 
gauge, A = t;(—Bx 2 , Bx 1: 0). We again use the decomposition into transverse 
modes, 

oo 

Hi = 0^>4, (4.1) 



n=l 



d 1„\ 2 { . d l,-,\ 2 fixn\ 2 



The first two terms at the r.h.s. denoted as h B describe a two-dimensional 
particle in the perpendicular homogeneous field. The resolvent kernel of such 
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an operator is well known ||14|| : 

(h B -zy\x,x') = 

x r 



1 (%B , , j \B\, 

— exp — (— xi% 2 + x 2 x x ) \x — x 

An V 2 ~ 4 



l\2 



\B\ 



z 



2\B\ 



U 



B\-z 



1; 



\B\ 



2\B\ ' ' 2 



\x — x 



where U is the irregular confluent hypergeometric function || 13.1.33]. For 
the sake of brevity we denote the exponential term in the above formula as 
& B (x, x'). The decomposition (|4.1|) then yields the sought resolvent kernel 



(H B -z) 1 (x, y; x', y') 



1 fiB i i , / \ \B\ I /|2 

— - exp I — {-x 1 x 2 + x 2 x 1 ) - — \x - x | 



n=l 
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C7 



\B\ - k 2 n (z) \B[ 

2\B\ ' ' 2 {X X 



'|2 



'nny\ / nny' 
x sin I — — sin 



4.2 The perturbed resolvent 



(4.2) 



As in the non-magnetic case we start from a single point interaction located 
at the point a G S and modify to the present situation the argument of 
Sec. [2.3| . We employ the fact that locally the magnetic field means a regular 
perturbation of the Schrodinger equation; motivated by this we define the 
one-center Hamiltonian in the same way as in Chap. 2: it acts as free 



(H*(a,a)iP)(x) 



(4.3) 



d 1 V / 9 



l 



<9 2 



9xi ' 2 ' ) \ dx 2 2 1 J dy 2 
for x a with the domain changed to 

D(H B (a,a)) = { ip G L 2 (£ \ {a}) : H B (a, a)ip G L 2 (£ \ {a}), 
-0(x, 0) = ?/>(x, c?) = 0, Li(ip, a) — 47raL ( , ?/ ; , a) 



(4.4) 



0} 



where Lq and Li in the last relation are again the generalized boundary 
values from || Chap. 1.1]. 



31 



Remark 4.1 To justify such a definition we have to check that the resolvent 
kernel (|4.2|) has the same singularity as the non-magnetic expression (|2.3|) for 



x' — ► x. We have Gq(x, x'; z) ~ c|af — by fl5| , Thm. III. 5.1] with a non- 
zero c independent of z. To check that the constant has the needed value it 
is sufficient to find operators H- and Hp with (3 from some interval (/?o, 1) 
whose resolvent kernels are ~ (47r) _1 |x / — x\ around the singularity and which 
satisfy the inequalities 

H-<H* </T%, (4.5) 
since the last named property implies easily 

for a fixed z < and c = (47r) _1 follows by contradiction. For the lower bound 
we choose the projection to the layer of the magnetic Schrodinger operator 
(-zV - A) 2 in L 2 (1R 3 ) obtained by removing the Dirichlet boundaries of E, 
because the latter is of the form hs ® I + I ® (—d 2 ) and we may apply 
the bracketing argument [[7], Sec. XIII. 15] to the non-magnetic part in the 



y-direction. Its kernel is known |16|, [17| to be 

G*(x,x';z) = ±fi|lY Vfos') ( 4 -6) 



v exp [-(\B\(2l + l)-zfl 2 \y-y'\] (B 

I ~t~ O Ol Dl 



/|2 

: ,„ J^/ I — IX — X 

\V 2 ( ^2' 

2 2|B| / 

where the factor $ B is the same as in the relation ft4.2| ) and are the 
Laguerre polynomials. We want to prove that \x — x'\Gq{ ) tends to 
(47r) _1 in the limit \x — x'\ — > 0. Putting |x — x'\ = 0, we can neglect 
the Laguerre polynomials and the factor $ B ; it remains to compute the 
simplified sum. For a given z there exists an integer number Z such that 
|-B|(2/ — 1) — z > 0. Then one could split the series into two parts: the 
finite sum over I = 0, . . . , lo — 1 which is irrelevant for the singularity and the 
truncated series with I = lo, . . . . The latter can be estimated as follows, 

L<M " a p[-(| fl |(2 1+ l)-^] < 

4 "6: (|bi(2/+i)- 2 ) i/2 - 

where 

\B\ f°° exp [-(|B|(2Z=f 1) -zfl 2 o\ 



dl = e 

An J, (\B\(2l T l)-zy/ 2 Attq 



-(|B|(2* =Fl)-^) 1/2 e 
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Hence the resolvent Gq has the needed singularity. 

In he opposite direction we add a Dirichlet boundary at \x\ = R cutting 
thus a finite cylinder of E. It is clearly sufficient to find a bound of the type 
(|4.5| ) for the "planar" part of the operator and this task is further reduced 
to finding bounds for its partial-wave components, 

~ B ld_ fm Br\ 2 

on L 2 (0,R) with the boundary condition lim r ^o0( r )( rl ' /2 l nr )~ 1 — at the 
origin |2], Sec. 1.5] and Dirichlet at r = R. A comparison of the potential 

/ \ 1/2 

terms shows that h 1 ^ < ft hP m holds if R < ( r|r (/3 -1 / 2 — 1)J , so one can 
choose for Hp the non-magnetic Dirichlet Laplacian in the cylinder of an 
arbitrarily small radius. This operator has again the resolvent kernel with 
the needed singularity ||. Notice finally that alternative ways to prove this 
result can be found in |18| or derived by techniques from the classical theory 



of partial differential equations [19, Thm. 20.6]. 



Remark 4.2 The previous remark still does not answer the question about 
the Green's function singularity fully. To explain why it is the case, recall 
that the requirement of symmetry of the magnetic Hamiltonian yields for 
any functions ipi, ip2 from the domain by means of the Gauss theorem the 
condition 




- dip 2 d^i - Ar 

--01— VWi~^ h2z^i^ 2 — I db 

or or r 



where f = x — a, r = \r\ and the integral is taken over the surface of the 
sphere Sg with center at a and radius r. It is satisfied if the functions have 
the following asymptotic behaviour in the vicinity of the point interaction: 



ip(x) 



1 + A(a)(x — a) 







\x - 


- a\ 



+ ci + 0(\x-a\). 



(4.9) 



This motivates us to change the generalized boundary value L\ to 

1 + iA{a){x — a) 



Li(ij},a) 



lim 

\3-a\-*0 



i/j(x) - L Q (ij},a) 



\x — a 



(4.10) 
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This suggests to use for £(£?, a; z) the following limit 



lim 

x— >a 



(£, a; z) 



1 + iA(a)(x — a) 
4tt\x — a| 



(4.11] 



which has the disadvantage that it is direction-dependent and can be altered 
by a gauge change. 

However, it is possible to employ the function £(-B, a; z) defined by means 



of the pole singularity alone. By Q4.2|) the Green function G B (x, a; z) has the 



form exp(iA(a)(x — a))F(\x — a|; z) in the symmetric gauge. In Remark |4J 
we have shown that F(\x — a\; z) has the following asymptotic behaviour for 
small \x — a|, 



F(\x 



1 



a\; z 



4tt\x 



+ c{z) + 0{\x-a\), 



where c(z) 



lim 



The function c defined in this way can 

be understood as £(B,a;z) obtained for the Green's function without the 
exponential factor; the corresponding generalized boundary value L\ does 
not contain the extra part A(a)(x — a). In the respective asymptotic formula 
for complete Green's function one has to multiply the last expression by two 
first terms of Taylor series of the exponential factor, 



Gn (x, a: z) 



1 + iA(a)(x — a) 
4tt\x — a\ 



+ c(z)+0{\x 



Using now the modified boundary value L x we find that £,(B,a;z) = c(z). 
This justifies the choice of the generalized boundary values made in the be- 
ginning of this section; one has to keep in mind to ignore the exponential 
term a) when computing £(-B, a; z). This argument remains also valid 

for a finite number of point interactions, hence the functions £(B,aj;z) for 
j = 1, . . . , N contained in matrix A(B, a, a; z) in Section |4.4| below are com- 
puted in the way described here. 

After this digression we can return to evaluation of the Green's function 
of the operator H B (a, a). By construction it is given by the Krein's formula, 



G B (x,x' 



G B i 



00 - *)C - z ^ I 



G B (x, a; z)G B (a } x'; z) 



€(B,a; z) 



(4.12) 
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where the regularized Green's function £(B,a; •) is, as we have explained in 
the above remarks, now given by 



£(B, a: z) = lim 
\s-s\-*o 



Gq (x, a; z) 



1 



An\x — a 



lim 



^ oo 



2nd 



n=i 



\B\-kl(z) 
2\B\ 

nnb 



U 



\B\ 



kl{z) \B\ 2 
2\B\ ' ' 2 8 



x sin 



d J 4:WQ 

One can check directly the consistency requirement 

G(x, x'\ z) 



lim G (x, x'\ z) 



(4.13) 



(4.14) 



for fixed a, a and x ^ x', where G(-, •; z) is the non-magnetic Green's function 
of Sec. p.3| . It follows easily from a known relation |5], 13.3.3] for the confluent 
hypergeometric function, 



lim T(u)U (u, 1,-)= 2K (2y/s) . 



To make use of the Green's function, we have to evaluate the r.h.s. of 
(|4~l3| ). We employ again the same trick and split a part of the series which 
can be summed explicitly for a general g while in the remaining series the 
limit g — > can be interchanged with the sum. Since U(u, 1; •) and Ko(-) 
have both a logarithmic singularity at zero, we modify the Ansatz of Sec. |2l| 
writing £,(B, a; z) = £i + £ 2 with 



im 

-*o 2nd ^ 

n=l 



\B 




- k 2 n (z) 


2\B 





u 



B\~kl{z) 
2\B\ 



1 \ B \ - 



nnb 
~cT 



6 



lim 



1 °° / 
_£2if (^)sin*( 

n=l x 



nnb 



2nd ^— f \ d J \ d j Ang 
The function £ 2 is evaluated as above: in analogy with (|2.13 ) we have 

6 =jl + r^C'V-ootC^v (4.i5) 
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As for the first part £1, we employ the small g asymptotics for the confluent 
hypergeometric and Macdonald functions, 



U{u,l;s) 
Ms) 



>s + ^(«) -2^(1)] +0(\slns\) 



as 







ln--^(l)j +0(, 2 ) 



as 



2T]| , 6.7.13] and || 9.6.13]). Putting them together we see that the 



[see 

summand behaves for small g as 

(\B\-z+(^ff 



In 



2|fi|d 2 



2151 



0(|^ 2 ln^|). 



For large n the digamma function i(j(s) = ln(s) — 2/s + 0(s 2 ), so the above 
expression can be written for large n as 



In l + {\B\-z) 



d 

irn 



2\B\ 



2 \B\-z+m 



+ 0(n- A ) 



d 

7m 



+ 0(n 



Hence the series in the definition of £i converges uniformly and the limit can 
be interchanged with the sum giving the sought formula for the regularized 
Green's function, 



2nd 



E 



71=1 



In 



7m 



+ 



And 



7 + t/> 



2\B\d 2 
b 



4 



\B\-z + {^y 



2\B\ 



sm 



d 



7T 

+ -COt 



lib 

~d 



irnb 
~d~ 

(4.16) 



Remark 4.3 The scaling behavior for the family S CT = IR 2 x [0, da], a > 0, 
is similar to that of Remark |2TI , however, one has to scale simultaneously 
the magnetic field by 

B a = a- 2 B. (4.17) 

In distinction to the previous sections we shall keep a general d in the fol- 
lowing discussion. 
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4.3 Spectral properties 

The essential spectrum of H B (a,a) remains the same as that of the free 
Hamiltonian H B which follows easily from Weyl's theorem |7|, Thm. XIII. 14]. 
The latter is in turn obtained from the essential spectrum of two-dimensional 
Landau Hamiltonian, u(h B ) = cr eS s(h B ) = { \B\(2m+ 1) : m G No} - see, 
e.g. Thm. 1]. Using the transverse- mode decomposition ( 4.1[ ) we arrive 
at 

a ess (H B ) = a ess (H B ) = ||fi|(2m + l)+ (™) 2 ■ m,n-lGN J. (4.18) 

Furthermore, the general properties of self-adjoint extensions mentioned in 
Sec. ^]4| imply that there is at most one eigenvalue in each spectral gap 
of the unperturbed operator, i.e. between the two neighboring values from 
&ess {H B ) and in the interval (—oo, \B\ + (n/d) 2 ); to find these eigenvalues 
one has to solve the equation 

t(B,a;z) = a. (4.19) 

We have already checked that the terms of the series expressing £(B,a;z) 
decay as n~ 2 so the series converges with the exceptions of the points where 
the ip function has a singularity, i.e. for any z G R \ a ess (H B ) . It is also 
obvious that £(£>, a; z) is real for any such z. We have 

>-> ^f>' ** (?) ■ <**> 

= 4-ip(s) is the trigamma function. For large n it behaves as 




B 


— z 




2\B 


2\B 




\B\-z + 



+ 0(n 



so the series converges for z G R \ cr ess (H B ) and we are allowed to inter- 
change the sum with the derivative. The explicit expression 6.4.10] for 
the digamma derivative, 

OO j 

^) = Er-^ for s ^ 0,-1,-2,... , 

p£ (s+J) 2 
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z 


( 





Figure 5: The function £(B, a; •) for £ = 1, d = 7t and b = it/6, n/3. 



shows that £(B, a; z) is monotonously growing with respect to z in each gap. 
Using the relation || 6.3.16], 

oo 

^(1 + a) = - 7 + V * s^-1,-2,..., 

we find that £ diverges as z approaches any point of the essential spectrum 
behaving in its vicinity as 

I Rl sin 2 ( 7rnb \ 
£(B,a;z) = - l -4 — 5 + 0(1) ( 4 - 21 ) 

provided the considered value from a ess (Ho) is "non-degenerate" in the sense 
that it can be expressed by means of a single pair of indices m, n. This is 
the generical case since the last name property is valid always if the ratio of 
the coefficients \B\ and (vr/<i) 2 is irrational. If it is rational, then to a given 
z G a ess (H ) there may exist different pairs mj,rij with the index belonging 
to a family J(z ) such that \B\(2rrij + 1) + (Trrij/d) 2 = Zq for all j G J(zq). 
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Taking this degeneracy into account we have 

\B\ 1 9 ( 7rn,-6\ „„. . 

« B - 3 ^ = -U— S sin hr) + ° (1 >- (422) 

When z approaches z from below, £ diverges to oo, while if it approaches 
zq from above, £ goes to — oo. The formula (|4.22 ) also shows that we can 



disregard those j for which sm(irnjb/d) = 0: as above the system does not 
"feel" a point perturbation situated at a transverse eigenfunction node. 

To find the behaviour as z — > — oo, recall that the argument leading to 
the expression ( |4. 16|) shows that the latter differs from the non-magnetic 
formula ( 2.14] ) by the replacement z — > z — \B\ together with the addition of 



terms which remain bounded as z — > — oo. Consequently, the asymptotics is 
independent of B and given by the formula 



£(B,a;z) = -- + 0(1) as z -»• -oo. (4.23) 

An 

Having discussed the properties of the function £, we can apply the con- 
clusions on the equation ( 4.19|) . Since £ is strictly increasing between every 



pair of neighbouring singularities, there is a unique eigenvalue e r (B,a,a) in 
each gap for any aeK and it satisfies the inequality 

e r (B,a,a) > e r (B,a',a) if a>oc, (4.24) 

where the index r is labeling the gaps of Hq. The value r = corresponds 
to the interval (— oo, \B\ + (n/d) 2 ), r = 1 represents the first finite gap with 
the left endpoint \B\ + (n/d) 2 , etc. 

Since we know the behaviour of £ around the singularities, we can write 
the explicit expression of s r (B, a, a) in the limits of strong and weak coupling. 
Suppose that the point \B\(2rrij + 1) + (nrij/d) 2 = z e a ess (H ) separates 
the r-th and (r + l)-th gap, then using the implicit-function theorem we get 

to -"\ 1 \ B \ sr- ■ 2 ( Krijb 
£ r (B,a,a) = z sin I — -j- 

+ £>(oT 2 ) as a^oo, (4.25) 



, ~ 1 IJBI v-^ 2 f KTljb 

e r+1 (B,a,a) = ab - ^Bm^— 

ieJ(zo) v 
+ 0(a~ 2 ) as a — > — oo 
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We see that the strong and weak limits are similar in the magnetic case. 
A different behaviour we find only for the lowest eigenvalue for which the 
asymptotic formula Q4.23 ) gives 

e (B,a,a) = -16vrV(l + 0(a -1 )) as a -> -oo . (4.26) 

A finer estimate with the error term replaced by 0(a~ 2 ) can be obtained 
when the £ is expressed in terms of the Hurwitz £- function - see, e.g., (17 . 
As in the non-magnetic case, the residuum of the resolvent pole in the formula 
( |4.12| ) yields a (non-normalized) eigenfunction 

ip(x; a, a) = G$ (x, a; e r (B, a, a)) (4.27) 

corresponding to e r (B, a, a), where Gq is the free resolvent (|4.2|). 

We are naturally interested in the behaviour of the eigenfunction in the 
limit \a\ — > oo for eigenvalues satisfying the asymptotic relations ( [4.25 ). 
In the both cases, by |2(], 1.17.11] the gamma functions in the sum ( f4.2| ) 
corresponding to rij, j G J(zq) diverges as 

(-l) m 3 1 

r(-m,- + e) = - — +0(1) as e^O. 

rrijl e 

This makes it possible to write the leading term of the asymptotics explicitly, 

ip(x;a,a) = 




d 

By H 13.6.27] the hypergeometric function is reduced to a Laguerre polyno- 
mial, 

U(-n,l,u) = (-l) n n\L n (u) , 

at positive-integer values of the first index, so the wave function can be 
rewritten as 

ip(x;a,a) w — - — —® B (x,a) (4.29) 

E ie j(, o) sin 2 (^ 



Lm J 

jGJ(^o) 



\B\ 2 \ / 7mjb\ (Ttnjy 

\x — a\ sin — - — sin ' 



2 1 1 J \ d J V d 
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In the last formula we have not specified the error term which has an ex- 
tra part coming from the variation of the hypergeometric function U(-, l,u) 
around integer values. 

We stated that the eigenvalue Eq(B, a, a) in the strong coupling limit case 
has the same behavior as if there were no magnetic field. The same is true 
for the corresponding eigenfunction which is strongly localized, 

if>(x;a,a) = ~ a) ^ Un^j \B\ + 16n 2 a 2 + (^) V - a\ \ (4-30) 

Let us summarize the results obtained for the layer in the magnetic field 
with one point interaction. 

Theorem 4.4 Let H B (a,a) be defined by (\4-3j) and (\j-4\ ); then 

(a) a ess (H B (a,a)) = {\B\(2m + 1) + (rm/d) 2 : m,n-lGN }. 

(b) For any a£l there exists a single eigenvalue e r (B,a,a) between every 
two neighbouring values from the essential spectrum, with the exception of the 
case when the leading term coefficient in (U.2Q ) is zero and the eigenvalue 
coincides with the Landau in question for a particular value of a. Each of 
these eigenvalues is increasing and infinitely differentiate w.r.t. a. The 
corresponding eigenfunction is given by ([Jl 



(c) In the both limits a — > ±oo the eigenvalue behaves similarly according to 
( \4-.23j ), except for the strong coupling limit for the lowest eigenvalue, where the 
formula (\j-2b\ ) is valid. The eigenfunctions corresponding to the eigenvalues 



( \4-23j ) and ( \j-2dj ) are given by ( \j-28j ) and ( U-3Q ), respectively. 



4.4 Finite number of point interactions 

Next we consider a finite number iV of point interactions supported by points 
Sj, j = 1,...,N. We define the Hamiltonian in a way similar to that of 
Sec. 3, i.e. it will be given by ( fO| ) and (^4.4j), where a, a are again shorthands 
for a = (d\, . . . , on) and a = (a±, . . . , ajy), and instead of a single boundary 
condition there is an iV-tuple of them, 

L 1 {i/j } a j )-4nra j L {i/j ) S j ) = } j = l,...,N. (4.31) 



41 



Accordingly, the Krein's formula reads 
(H B (a,a) - z)~ x {x,x*) 



N 



= Gq(x,x'; z) + XjkjB, a, a; z) Gq (x, ay z) Gq{cl^ x*;z). (4.32) 
j,k=i 

Repeating the argument of Sec. 3.2, we find X(B, a, a; z)~ l = A(B,a,a; z) 
with 

A mj (B, a, a; z) := 5 jm (a rn - £(B, a m ; z)) - (1 - 5 jm )Go(a m , ay, z) , (4.33) 
or more explicitly by means of (|4.16|) and ( |4.2|) 

\B\-k 2 n (z) 



A 



.13 



2nd 



n=l 



111 



\2\B\d?) 



And 



7 + ( Uj i ) + ~ cot ( ^ 



-1P 



nb q 



2\B\ 



sin 



(4.34) 



and 



A m j 



2nd 

x C/ 



1 °° / 

s * B («-«,)E r ( 

n=l v 



1^1 

2|5| 



(4.35) 



1*1 -*M . |5| 2 . 

, 1; — — \a m -ay I sin 



2 B 



nnbj 
~d~ 



sm 



for j 7^ m. If some perturbations are arranged vertically, aj = a m for j ^ m, 
the last expression can be written as 



A 



in j 



1 

~ 2n~d 
x sin 



In 



£ 

n=l 

nnb m 
d 



/m 2 \ .( \B\-z+(^y 

\2\B\d?) W \ 2\B\ 



(4.36) 



sm 



nnb. 



+ 6 



|6j - br, 



as we find by a direct modification of the argument yielding £(B, a; z) 
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Consider again a point z e a ess (H Q ). If z approaches this value, the 
appropriate contributions to the sums in (|4.16| ) and (|4.2|) diverge and the 
matrix elements of A behave in the limit z — > zq as 



M V sin 2 — p + O 1 , 4.37 



&j Q m | 



ieJ(z ) 

\B\ 1 

x E (zi£L g ^l + (?)'-» |1; jgl 

. f urijbj \ . f imib m \ 
x sin I — J sin I — - — J + 0(1) if j m, 

ieJ( 2 o) v / v / 
if jV™. Oj = a m . 

Moreover, the hypergeometric functions reduce to the Laguerre polynomials 
in the limit, so 

\B\ 1 

A mj « _l_J$B( amjaj .) (4.38) 

l-Bl . 2 \ . { imibj\ . ( TTTlibr, 

x > I — \a m — aj] I sin I — r-i I sm ' 



ieJ(zo) 



ti / V d 



For large negative energies, z — > — oo, we employ the asymptotic behavior of 
£ given by ([4.23 ) arriving thus at 



Aii - «i + ^ + ^(l) (4.39) 

with .B-dependence being hidden in the error term. The non-diagonal part 
of the matrix A vanishes in the limit of large negative energy, because by || 
13.3.3] one has 



^ oo 

Ami « - $B( " m ' ^ S ^° (V|£| - &n 

n=l 



7mbj\ . / 7rn6 r 
x sm I — sm 



d 
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which is exponentially decreasing, K (s) « \pf s &~ s as s — > oo . In the case 
of a vertical arrangement, = a m for j 7^ m, we can again claim only that 
the matrix element Aj m is bounded. 

Having obtained the coefficient matrix A(B, a, a; z) given by ( 4.34|) and 
(|4.35| ) we proceed with finding the eigenvalues (the essential spectrum is 
preserved, of course). In the same way as in Sec. [3.3| we check that an 
eigenvalue z is a solution of the implicit equation 

det A{B,a,a;z) = (4.40) 

and the corresponding eigenfunction can be written as 

TV 

= d i G o(X' %; z ) > ( 4 - 41 ) 

i=i 

where d := (di, . . . , d^) is an eigenvector of the matrix A(B, a, a; z). 

As in the Sec.^]^ we would like to say something about the number of 
eigenvalues in each gap. We already know that this number is limited from 
above by N. Denoting a.j = dtj + a we can split the matrix A(z, a) into two 
parts 

A(z,a) = aI-M(z,a), (4.42) 

where a := (c?i, . . . , 5jv). The explicit form of the matrix M can be obtained 
from the relations (|4.34j) and (|4.35| ) by changing the signs and substituting 
aij by a,j. The same can be done for the formulae ( 4.37|) which express the 
behavior around the singularity. In the limit z — > Zq one can neglect the 
parameters a/s, then it is possible to write 

M (z, a) = S M(z ) + 0(1) , 



zq — z nd 



where 



sin 

ieJ(zo) 



2 / uriibj 



d 



M mj (z ) = ^ B (d m ,d j ) ^2 L r 

ieJ(zo) 



\B\ 

T 



' Txriibj\ ( irriibr, 
x sm 1 — I sin 
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Following [^TJ the eigenvalues of A(z) monotonously increase between two 
neighbouring singularities. Then all eigenvalues of M are non-negative num- 
bers. For 52j =1 J2i<=j(z ) sin 2 (irriibj / d) ^ at least one is positive, hence at 
least one eigenvalue of the matrix M(z, a) diverge to — oo or to oo as z — > zo + , 
z — > z _ respectively. Let us summarize the results: 

Theorem 4.5 Let the operator H B (a, a) be defined by ( \4-3\ ) and ( U-A ); where 
a — (Si, ... , a N ) and a = (cti, . . . , a^) and the boundary condition is replaced 
byfitfffy; then 

(a) a ess (H B (a, a)) = {\B\(2m + 1) + (vm/d) 2 : m,n-lGN }. 

(b) For any a G M. N there exists at most N eigenvaluee between every two 
neighbouring values from the essential spectrum with the multiplicity taken 
into account. The corresponding eigenf unctions are given by (\j.4- A) - 

(c) In the limits mini<j<7vaj — > oo and maxi<j<7vaj — > —oo at least one 
eigenvalue converges to each value zq from the essential spectrum from below 
and from above, respectively, with the exception of the case when the leading 
term coefficients in ( \j.3Tj ) are zero. In the strong limit case there are also 
N eigenvalues given by j\4 . 26j ) with corresponding eigenf unctions given by 
( \j.3C\) , where a and a is replaced by otj and Sj for j = 1 . . . N. 

Further results on the number of eigenvalues of the Hamiltonian Hq can be 



obtained in the same way as in p3|. 



Remark 4.6 If some of the point interactions are vertically arranged, we 
cannot exclude that eigenvalues are absent in a particular gap for some a. 

Consider two point interactions placed at (0, 0, l/4d) and (0, 0, 2/3d) with 
ol\ = 0L2 = 0. The numerical calculation for B = 1 and d = tt shows that 
the eigenvalues of matrix M(z,a) for z G (;z (l, 1), -2o(0, 2)) cover whole R 
except one gap, see Fig. |5|. The symbol z (m, n) represents the Landau level 
z = \B\(2m + 1) + (rcn/d) 2 . Hence for a from this gap the matrix A has no 
eigenvalue in the interval (zq(1, 1), zq(0, 2)) = (4, 5). 



5 Conclusions 

We have analyzed here spectral and scattering properties of a hard-wall layer 
with a finite number of point interactions. The results offer one more illus- 
tration of efficiency of Krein's formula which allows to reduce the task in 
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Figure 6: The behavior of two eigenvalues of the matrix M(z, a) for z G (4, 5). 



fact to an algebraic problem. There are other interesting question related to 
systems with finitely many perturbations such as relations between the per- 
turbation configurations and the spectra including the nodal structure, etc., 
positions of resonances including those coming from perturbation of the em- 
bedded eigenvalues, and so on. To keep this paper within reasonable limits, 
however, we postpone these questions to a sequel. 

The same applies to systems with an infinite number of point obstacles 
which offer a wider variety of spectral types. Let us briefly mention several 
problems which we regard as worth of attention. One of them concerns the 
number of gaps in periodic systems. A periodic layer of point interactions 
in M 3 has at most one gap ||, Sec. III. 1] . On the other hand it is known 
that the presence of boundaries can enhance the number of gaps in the two- 
dimensional case ||; a similar effect is expected in dimension three: for a 
thin enough layer there will be many open gaps between the first and the 
second transverse threshold. A more difficult question concerns the validity 
of the Bohr-Sommerfeld conjecture in such systems. 

Even more interesting are spectral properties of periodically perturbed 
layers in presence of a magnetic field. It is well known that a combination 
of a square lattice of point interactions and a homogeneous magnetic field 
leads to a very rich spectrum whose properties depend substantially on the 
number-theoretical properties of the ratio between the lattice spacing and the 
field intensity (which determines the cyclotronic radius) - see @, Sec. III. 2. 5] 
or pi ]. Putting such a system into a layer brings a third parameter (the layer 
width d) which will determine how "thickly" the transverse-mode component 
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are overlayed in the spectrum. 

The same applies to edge-type states. It was shown recently that an 
equidistant array of point interaction in combination with a homogeneous 
magnetic field can produce bands of absolutely continuous spectrum away 
of the Landau levels f25fl . One is naturally interested how the spectrum will 
change if the array is confined between a pair of hard walls. Other problems 
concern aperiodic perturbations, external electric field, spin effects, etc. 
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